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ABSTRACT 



This thesis is a report of an experimental investigation of the 
influence of horizontal mechanically induced simple harmonic vibrations 
upon the rate of heat transfer by free convection from a horizontal heated 
cylinder to air. The diameter of the cylinder was 7 / 8 "^ the ranges of 
the principal experimental variables were as follows: 



temperature potential^ 
frequency of vibration^ 
amplitude of vibration^ 
intensity of vibration^ 



d:.t: 150 °F ± 1 °F 

f; 104 cps 
a: 0.025 to 0.128" 

af: 0.21 to 1.11 ft/sec 



The data show that the influence of vibrations upon the heat- 
transfer coefficient is negligible for intensities of vibration^ af^ less 
than 0.3 ft/sec; above this so-called "critical intensity^" the effect of 
vibrations is to increase the heat-transfer coefficient significantly. 

The quantitative heat transfer results of this investigation show close 
agreement with the results reported for horizontal cylinders in the pre- 
sence of high intensity horizontal sound fields . An elementary flow- 
visualization examination utilizing smoke was performed which indicated 
that the free convective flow near the lower test cylinder stagnation 
point remained laminar in the presence of vibrations. The quantitative heat 
transfer results, supported by the flow study, lead to the conclusion that 
the physical mechanism of heat transfer for mechanically induced horizontal 
vibration of a cylinder is the same as for acoustically induced horizontal 
vibrations in the fluid medium surrounding the cylinder. The close quanti- 
tative agreement of the results for the two methods of vibration also indi- 
cates that the influence of the amplitude of vibration to cylinder diameter 
ratio, u/Dq, which differed by a factor of 12 . 3 ^ is not significant for the 
range of a/D^ ratios considered. 

Thesis Supervisor: George A. Brown 

Title: Assistant Professor of Mechanical Engineering 
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CHAPTER I 
INTRODUCTION 



The study of heat transfer in the presence of vibrations is a rela- 
tively new area of investigation that has received increasing attention in 
recent years. Several experimental investigations have indicated a marked 
increase in convective heat-transfer rates in the presence of strong vibra- 
tions . Since heat transfer is a limiting factor in many areas of engineer- 
ing design, it is important to account for the effects of vibration on 
heat transfer where vibration is present. In certain cases, advantage may 
be taken of increased heat-transfer rates to render certain heat-transfer 
dependent processes more efficient. As an example, Raben and Dietert 
have conducted an investigation designed to use vibrations to increase 
heat-transfer rates in a saline water conversion process. 

The influence of vibrations upon the free convective heat-transfer 
rate from a heated surface can be investigated experimentally by two dif- 
ferent methods. In one method the heated surface is held stationary 
within a fluid medium. An oscillatory motion is impressed upon this fluid 
medium -- that is, this method requires the establishment of a sound field 

1 

in the vicinity of the heated surface and will, therefore, be referred to 

I 

subsequently as the "sound method." An alternate experimental approach is 
to impress a vibratory motion on the heated surface itself, leaving the fluid 
medium surrounding the surface otherwise undisturbed. This second method 
generally employs a resonant mechanical system to achieve sustained vibra- 
tions of the heat-transfer surface and will, therefore, be referred to as 
the "mechanical vibrations method." The objective achieved by the two 
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methods is identical. That is^ the relative velocity vector "between the 

fluid medium and the heated surface undergoes a periodic change of magnitude 

(o) 

and sign. Indeed^ an analytical study^ ^ has shown that^ for sufficiently 
long wave-lengths^ identical fluid flows should exist near a surface if the 
direction^ amplitude, and frequency of vi"bration are the same, regardless of 
which of the two methods is used to create the vibration. 

An early attempt to correlate the effects of vibration on heat trans- 
fer is that of Martinelli and Boelter^^^. In this investigation, a SA" 
diameter horizontal electrically heated tube, immersed in water^ was sub- 
jected to vertical sinusoidal vibration. Amplitudes up to 0.10" and fre- 
quencies up to 40 cps were employed. The range of temperature potential, 

^t, extended from 8 to 45^F. No influence on the heat-transfer coefficient 
was observed in the lower range of vibration Reynolds number; this result 
was attributed to the dominance of free convection in this range. At higher 
vibration Reynolds numbers, however, a gain' of as much as 400 per cent rela- 
tive to the non-vibration case was observed. Unfortunately, Boelter, in a 

( 4 ) 

communication to Lemlich , has reported that these results could not be 
verified in a later series of experiments . 

The effect of vibrations on convective heat-transfer rates from 
horizontal wires to air has been studied by Lemlich^^^. Three heated 
wires of varying diameter (0.0253”, O.O396" and O.O8IO") were set in vibra- 
tion at frequencies between 39 and 122 cps with a maximum amplitude of 0.115”* 
Increases in the heat-transfer rates were reported for increases in frequency 
or amplitude. Improvements in heat-transfer coefficients by as much as a 
factor of four were reported. In addition, Lemlich reported similar results 
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for both vertical and horizontal vibrations and concluded that the varia- 
tion in heat transfer was independent of the direction of vibration. 

Tsui^^^ has conducted an experimental investigation of the effect 
of transverse vibrations on the heat-transfer rate from a heated vertical 
plane surface. Increases in heat-transfer coefficients up to 25 per cent 
were noted. Tsui attempted to treat the system analytically, but conclud- 
ed that the general equations correlating vibration and heat transfer are 
insolvable with present techniques. Shine^^^ performed experiments similar 
to those of Tsui. A heated vertical plate was subjected to transverse vi- 
brations of 11 to 315 cps at amplitudes up to O.061". Plate temperature 
range was from I3I to 2 T 9 °F* Significantly, these tests demonstrated that 
as the amplitude -frequency product, af, (arbitrarily defined as the vibra- 
tion intensity) was increased, the heat transfer was not affected until a 
specific value of af was reached. Above this "critical" vibration inten-. 
sity, af = 0.083 ft/sec, the rate of heat transfer increased markedly. 

(t) 

Freedman ' has investigated the case of a vertical heated cylinder 
vibrating in an axial direction. No change in heat-transfer rate was 
noted over the range of vibration intensities achieved (maximum vibration 
intensity, af = 0.2 ft/sec). 

Fand, Kaye and Teleki^®-’ have recently conducted experimental 
investigations of the influence of vertical mechanical vibrations 
(Fig. 1 -b) upon the heat-transfer rate by free convection from a heated 
horizontal cylinder to air. The diameter of the cylinder was T/ 8 ", and 
the ranges of the variables were as follows: temperature potential, At, 

25 to 185 °F; amplitude of vibration, a, 0 to O.16"; frequency of vibration. 
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54 to 225 cps; and intensity of vibration^ af^ 0 to 1.22 ft/sec. These 
authors report that the sole controlling vibrational variable is the vibra- 
tion intensity^ af. Further^ the critical value of af was about O .3 ft/sec. 
A flow- visualization study using smoke indicated that the fluid-dynamical 
mechanism which causes the increase in heat transfer rate is vibrationally 
induced turbulence . 

An extensive investigation involving the sound method^ that is^ the 
effect of sound fields upon the heat transfer by free convection from a 
stationary horizontal cylinder to air^ has been reported by Fand and 
Kaye^^^^. In their experimental arrangement the sound vibration vector 
was normal to both the cylinder axis and the gravity field (Fig. 1-a). 

The frequency range covered was from 1100 to 6l20 cps; the vibration in- 
tensity range was from 0 to I .3 ft/sec; and the temperature potential was 
varied from 0 to 250^F. A critical vibration intensity of approximately 
0.3 ft/sec was observed; below this intensity the heat-transfer coefficient 
was unaffected. The data also showed that for sound waves whose half wave- 
length was six or more times the test cylinder diameter^ the heat-transfer 
coefficient is a function of only two variables; namely^ the temperature 
potential^ >4t^ and the vibration intensity^ af. In the tests above the 
critical intensity the rate of heat transfer increased rapidly until an 
af of about 0.6 ft/sec was reached; above this intensity level the rate of 
change of heat transfer was less rapid and could be expressed by the fol- 
lowing algebraic expression: 
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0.722 



A t (af)^ 



1/3 



( 1 ) 



where is the coefficient of heat transfer in the presence of vibra- 
tions^ Btu/hr-ft^-^F; af is the vibration intensity^ ft/sec; and /it 

o 

is the temperature potential^ F. 

Fand and Kaye^^^^ examined the physical mechanism which caused the 
heat-transfer rate to increase in the presence of sound by means of a 
smoke flow- visualization technique. This study showed that at the criti- 
cal sound intensity the typical free -convective boiondary-layer flow pat- 
tern aroTind the heated cylinder is disrupted^ and a fundamentally differ- 
ent type of bo\indary-layer flow^ called thermoacoustic streaming^develops . 
This thermoacoustic streaming appears as a special type of vortex motion. 
In the presence of both horizontal vibration and a temperature gradient^ 
a pair of vortices alternately form in the two quadrants above the heated 
cylinder surface. The inception of this oscillating vortex pair forma- 
tion occurs at the critical vibration intensity^ and the further develop- 
ment of these vortices corresponds to that range of af wherein the heat- 
transfer coefficient increases most rapidly. A further increase in the 
vibration intensity serves to more fully develop these vortices^ but does 
not alter their form. 

The investigations summarized above have^ for the most part, been 
of an experimental nature . All attempts to obtain a general analytical 
solution to the problem have been unsuccessful thus far. In view of this, 
it appears necessary to add to the knowledge of the physical relationships 
between vibrations and heat transfer by empirical means . The primary 
consideration involved in the experimental program "undertaken in this 
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work is to establish the heat-transfer behavior of horizontal cylinders 
under horizontal mechanically induced sinusoidal vibration, where the 
vibration vector and the force of gravity are in the same relationship to 
a heated test cylinder as in the sound investigations of Fand and Kaye • 

No direct comparison of the sound method and the mechanical vibrations 
method with both the test cylinder axis and the vibration vector normal 
to the direction of the gravity field (Fig. 1-a) has been reported. As 
noted above, investigations made of the influence of vertical vibrations 
on the heat-transfer rate of a horizontal cylinder to air (Fig. 1-b) have 
indicated that the basic mechanism of convective air flow around the cyl- 
inder is vibrationally induced turbulence^ not thermoacoustic streaming 
as observed in the sound method. Therefore, it was desired to determine 
if the character of heat transfer and the associated thermo acoustic stream 
ing phenomenon exhibited in the sound method can be related to the mechani 
cal vibrations method. Accordingly, this study is an experimental inves- 
tigation of the influence of horizontal mechanically induced simple har- 
monic vibrations upon the heat-transfer rate by free convection from a 
heated horizontal cylinder to air. 



I 
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CHAPTER II 

GENERAL DISCUSSION OF EXPERIMENTAL DESIGN 

The primary objective of this investigation was to determine experi- 
mentally the influence of transverse horizontal mechanically induced simple 
harmonic vibrations upon the heat-transfer rate by free convection from a 
heated horizontal cylinder to air. This configuration (Fig. 1-a) was 
chosen in order to provide direct comparison with available data for sta- 
tionary horizontal cylinders in the presence of transverse; ‘horizontal 
sound fields. In particiLLar, the experimental program was designed to 
establish whether or not a quantitative comparison of the heat-transfer 
variation for the two cases at identical temperature potentials, ^t, 
and intensities of vibration, af, does in fact exist. To this end, then, 
the following experimental goals were established: 

a. To construct a mechanical vibratory system which, at resonance, could 
drive the test cylinder at vibrational intensities, af, up to 

1.50 ft/sec 

b. To provide a suitable means to measure vibration amplitude and 
frequency accurately 

c. To construct a heat transfer test cylinder which could operate 
under the imposed vibration at surface temperatures up to 350°F 

d. To provide a power supply to heat the test cylinder and a tem- 
perature measiirlng system for accurate temperature measurement 



of the test cylinder surface . 
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Based on the experience of others it was determined that the 

test cylinder should be rigidly coupled to a resonant mechanical system. 

At resonance the force magnification factor could be utilized for obtain- 
ing maximum vibrational amplitude. Also^ although the test program con- 
templated here involved only horizontal vibrations^ it was considered 
economical to design the vibration test equipment with a view toward 
future work employing other vibration directions. Hence^ it was decided 
to incorporate into the vibration system means by which the vibration 
direction could be varied from 0 to 90 ^ with respect to the direction of 
the force of gravity. 

The basic element of the mechanical system chosen to provide vibra- 
tion was a simple free-free beam as shown in Fig. 2. The free-free beam 
offers the advantage of confining the major forces^ which are the inertia 
forces^ to the beam itself and not to. the supporting structure. Support- 
ing the free-free beam at the nodal points minimizes any hysteresis losses 
in the supports since there is no force transmission to the supports. 
Therefore^ less energy absorption occurs at the nodal supports than would 
occur; for example; at the supports of a clamp-clamp or hinge-hinge beam. 

The test cylinder was rigidly connected to the test cylinder by 
means of a yoke. The beam itself was supported at the nodes by a "girabal” 
arrangement attached to a rigid "alignment platform" as shown in Figs. 2 
and 3* The distance between the gimbal supports could be varied to the 
computed nodal positions of the beam when the beam was placed in resonant 
vibration. The alignment platform was attached to a MB Electronics Model 
GIF Vibration Exciter which provided the external driving force. The 
motion of the test cylinder resulting from an electric excitation of the 
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coupled dynamic system was a uniform, one -dimensional sinusoidal trans*- 
verse vibration. The frequency of vibration employed was 10^ cps; however, 
beams designed for other resonant frequencies could be used, as provision 
was made for fitting beams of variable dimensions into the vibration test 
stand. Vibration amplitude was measured by a cathe tome ter • 

The test cylinder, which was heated by an internal electrical heater, 
had a diameter of 7/8"* Test cylinders previously used in other experi- 

( 8 , 9 ) 

ments involving heat transfer under conditions of vibration en^loyed 

guard heaters on either side of the central test section heater in order 
to eliminate the problem of axial heat conduction away from the test sec- 
tion. The complexity resulting from this type of construction — that 
is, the use of three interconnecting parts, each with its own heating 
element and thermocouple arrangement -- has resulted in premature failure 
by fatigue. In order to avoid fatigue failure, the test cylinder for this 
project was designed with one heater; a metering device was installed at 
each end of the cylinder to measure the axial conduction heat losses to 
the test cylinder supports . 

A total of eight thermocouples were installed in the test cylinder 
for surface temperature measurement. The location of these thermocouples 
is as shown in Fig. 5* Each copper-constant an thermocouple was encapsu- 
lated in the end of a copper-plated stainless steel tube. The tube was 
then closely fitted into its individual hole drilled axially into the 
cylinder wall. Thus the thermocouples were firmly fixed in position, 
and a good thermal contact between the jiinctions and the cylinder was 



established. 
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There vas an obvious danger of test cylinder lead-wire breakage 
under vibrationally induced loads as high as 200 g's during the course 
of the tests. Several special techniques were used to avoid wire break- 
age and are described in detail in Chapter V, Section 1-c . 

The vibration test stand itself was located in a special test room 
8 X 13 X 15 ' high. A 4 X 6' glass window was placed in the wall between 
the test stand and the adjacent control room. Thus, optical amplitude 
measurements and all other test operations could be conducted without 
disturbing the air in the vicinity of the test cylinder. 
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CHAPTER III 
APPARATUS 

The apparatus can he conveniently divided into two basic categories 
with subdivisions as follows: 

1 . Vib r at io n Apparatus 



l-a. 


General 


1-t). 


Resonant beam 


1-c. 


Test stand 


1-d. 


Vibration exciter and exciter power supply 


1-e. 


Amplitude and frequency measuring apparatus 



2. Heat Transfer Apparatus 



2-a. 


General 


2-b. 


Test cylinder 


2-c. 


Power supply 


2-d. 


Temperature measuring apparatus 



1. Vibration Apparatus 
1-a. General 

In order to produce the Intensities of vibration required for this 
investigation^ it was necessary to employ a mechanical resonant system. 
This was evident from the characteristics of the available vibration 
exciter; the maximum acceleration obtainable from this exciter^ using 
direct force take-off^ was about 75 As stafced in Chapter II^ ac- 

celeration loads up to approximately 200 g’s were required to obtain the 
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desired value of af = 1.5* To oltain these high accelerations^ it was 
necessary to interpose a mechanical resonant system between the vibration 
exciter and the test cylinder. This is the general technique described 
in the literature used to produce high values of vibration intensity. 

An additional requirement of the vibration system was that the vi- 
bration motion produced had to be sinusoidal and one -dimensional in 
character. This requirement restricted the selection of the resonant 
system to one which would vibrate essentially in its first natural mode 
with negligible effect from any higher harmonics. A prismatical beam was 
the system chosen which satisfied both the resonance and sinusoidal^ one- 
dimensional motion requirements. 

Other mechanical resonant systems could have been devised to pro- 
duce the motion required — such as a system using heavy coiled springs^ 
for example. However^ the simple geometry of the rectangular beam lends 
itself more conveniently to accurate analysis^ design and construction. 

Resonaint beam vibration with clamp-clamp^ hinge-hinge and free -free 
boundary conditions were investigated. The free-free vibration of a beam 
in the first mode was selected as the preferable method for obtaining the 
required vibration intensity because this system exhibits the least damp- 
ing. In free-free beam vibration the major forces^ which are the inertia 
forces^ are balamced within the beam itself and not transmitted to the 
structure supporting the beam. Thus^ the joint type of damping loss is 
eliminated. An additional advantage results from this absence of major 
force transmission at the supports: the supporting structure is largely 

free of induced vibrations and can be designed accordingly. 
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One disadvantage in employing a free -free resonant beam compared to 
a clamp-clamp resonant beam is the necessity to provide a different beam 
for each vibration frequency desired. However, since it has been demon- 
strated by other investigators^®-’ that the sole vibrational variable 
affecting heat-transfer rates is the intensity of vibration — the ampli- 
tude-frequency product, af. Hence, it is possible to investigate the 
Influence of vibration intensity on heat transfer by varying the amplitude 
of vibration at constant frequency. 



1-b. Resonant beam 

The complete analysis of the free-free beam with center-span test 
cylinder load is given in detail in Appendix A. The results and interpre- 
tation of this analysis are presented here. 

In the analysis it was assumed that the beam at its resonant fre- 
quency, has a motion essentially of the first mode only, the effects of 
the higher harmonics being negligible. This assumption is considered 
valid for a vinlform, carefully machined, prismatical beam. Further, tests 
showed that the vibratory motion of the beam was, in fact, slmrsoidal and 
one -dimensional in character. 

The resonant frequency equation of the system, as derived in Appen- 
dix A, in terms of beam mass, mj^, and the mass at beam center-span, m^,, 
is: 

(sechof>^ + cos c</) 

(sinflf/’i- cos tan ho(/ ) 



yf is one-half the beam length and is defined as : 



k 



0^ 



12 X" 
Eg© 




where 



f 













► 



2 

E = modulus of elasticity, Ib/in 

2 

= gravitational constant, in/sec 
= beam weight density, Ib/in'^ 

tO = circular frequency, rad/sec 

h = beam height in motion direction, in 
Equation (2) is plotted in Fig. 6. 

In terms of the beam physical parameters, the resonant frequency 
can be expressed as: 



1 go 1/2 E 1/2 „ 



(3) 



Therefore, for a beam of given dimensions, the vibration intensity is: 



E 1/2 

af = Constant x a x (r;) 

IT 



(M 



Consideration of maximuin tending stress leads to the following rela- 
tionship between maximum bending stress in the beam^ ^max^ amplitude 
of vibration^ a^ as derived in Appendix A and plotted in Fig. 9* 



= I (a function of 0(/ ) (5) 

Equations (2), (4) and (5) form the essential basis for resonant 
free-free beam selection. Fig. 9 clearly shows that a high value of is 
desired in order to minimize bending stress in the beam^ and Fig. 6 indicates 
that large values of c</ correspond to large values of the beam mass to 
center load mass ratio ^ in^/m^. The weight of the center load was known to 



be approximately 5*5 lbs. This weight included the weights of all the 
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moving parts attached to the center of the heam: the test cylinder yoke^ 

the moving element of the vibration exciter^ and the parts used to couple 
the vibrator moving element and test cylinder yoke to the beam. There- 
fore^ it is clear from Fig. 6 that for a reasonable beajn size^ the value 
of woiild have to he within the range of about 2.20 to 2.28. These 



from Fig. 6, using a value for center-span load of 5*5 lbs. In the actual 



weight of 31*6 pounds were used for purposes of analysis. A lower value 



would require a heam larger than could he conveniently used. 

Equation (4) indicates that a high value of E/y^ is desired. 

Additional considerations involved in the heam design include thei 
vibration frequency and fatigue strength. The frequency at which reson- 
ance was desired was in the neighborhood of 100 cps. This frequency was 
selected as being high enough to he well within the frequency ranges used 
by other investigators to demonstrate marked effects of vibration on heat- 
transfer rates, and low enough to avoid more severe inertia loads at yet 
higher frequencies. The loading is proportional to the square of the 
frequency, and even at frequencies of 100 cps, inertia loads of about 
200 g's were anticipated. The vibration test stand was designed so that 
beams of varying dimensions could be Installed; hence, frequency could be 
varied by changing beam size. 

Fatigue strength was a primary factor involved in beam material 




from 19 lbs to 52 lbs as obtained 




2.25 and a corresponding beam 




would impose a greater bending stress, and a higher value of cX 



selection. The beam was required to sustain complete stress reversals 
for long periods of time . Beam failure might have resulted in test 
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cylinder destruction or, at least, would have caused considerable delay 
in the test program. 

Of the several beam materials considered, aluminum and steel were 
given consideration based on strength and fatigue characteristics, avail- 
ability, machinability, and the requirement of Equation ( 4 ), namely, a 
high ratio. From the stress standpoint, the best available aluminum 

alloy was 2024-T^, a special heat-treated alloy which has a yield stress 
of 50^000 Ib/in^, and the best available steel alloy was a special treat- 
ment steel, "HY-100," which has a yield stress in excess of 100,000 Ib/in . 

The endurance limit for A1 2024-T4 is quoted as 20,000 Ib/in^. This 
is an arbitrarily defined endurance limit, for it is well known that the 
fatigue characteristics of non-ferrous metals do not exhibit a limiting 
value of fatigue strength below which the material strength is independ- 
ent of the number of stress cycles. No fatigue data for HY-100 steel is 
currently available. It is known that the endurance limit for most 
alloy steels is approximately one-half the yield strength; therefore, 
’’HY-100" steel was considered to be adequate from the fatigue standpoint 
up to stress levels of approximately 50,000 Ib/in^. 

Each of these two metals could be machined readily; and aluminiom 

8 8 

and steel have about the same E/^ ratios -- 1.0 x 10 *^ and 1.1 x 10 
inches, respectively -- giving each material about the same footing rela- 
tive to the requirement of Equation ( 4 ). 

The width, b, of the beam was established as I. 50 " initially in order 
that existing test cylinder support equipment could be used. In a study 
made of A1 2024-T4 and HY-100 beams, the first comparison was made on the 
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basis of equivalent beam size^ and the second comparison was made for 
identical values of o</ . The results of the equal beam size calcula- 
tions made in accordance with the method given in Appendix B were: 



Beam Size -- 


1.50" X 


l.kk" X 


52,00'" 




of.i 


f(cps) 


af 


(ft/sec) Cr^ (ib/in^) 


steel 2.25 


100 




1.50 


52,000 


Aluminum 2 . 15 


100 




1.50 


25,500 


The results of the comparison 


on the 


basis of identical values of 


were : 










Size 




f(cps) 


af (ft/sec) 


^max 


steel 1.50"xl.i)-4"x52.00" 


2.25 


100 


1.50 


52,000 


Aluminum 1.50"x2.90"x73.00" 


2.25 


100 


1.50 


18,000 



These figures indicate^ for a/ = 2.25^ that either material would 

be satisfactory. Considering the superior fatigue characteristics of 
steel as compared to aluminum^ however^ HY-100 was selected as the beam 
material for the final design. 

The HY-100 obtained for the project had the following alloy content 
(MIL-S-2015i|.A): 

Carbon Nickel Chromium Molybdenum 



0.20-0.21 



3.50 



1.50 



0.40 
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Two standard tensile specimen tests gave the following results: 

No. 1 No. 2 

Yield strength (ih/in^) 106^000 102(^200 

Ultimate strength (ih/in^) 124,300 111,900 

The dimensions and characteristics of the beam used in the test 
stand were as follows: 

widths b height^ h lengthy 21 f (calc.) f (actual) 

1.4t4" 1.4o8" 50.00" 105.3 cps 104cps 

Small variations from the planned dimensions were necessary because of 
the machining required to produce a straight beam with a 63 micro-inch 
finish from the available stock. In addition^ the length of the beam was 
adjusted to raise the calculated resonant frequency somewhat above 100 cps 
to allow for the added mass of the gimbals mounted at the beam nodes and 
because the weight placed at beam center span was not accurately known 
imtil all parts were constructed. The measured resonant frequency was 
about 1 per cent less than the calculated value. 

1-c . Test stand 

The vibration test stand and its components are shown in Figs. 2, 

3 and 4. The purpose of the test stand was to provide a means of support 
for the resonant beam and a method for accurate alignment of the resonant 
beam with respect to the vibration exciter. 

Resonant beam support at the nodes of vibration was accomplished by 
two gimbal assemblies. At each node an aluminum gimbal frame was fitted 
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to the resonant heam by means of a stainless steel slipper. The beam 
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slipper was incorporated into the gimbal design so that beams of different 
sizes could be fitted into the gimbal frame — that is^ each beam with its 
own pair of slippers could rest in the existing gimbal frames. A beam 
clamp plate served to hold the beam firmly within the gimbal frame by means 
of two l/4” set screws and lock-nuts. To allow for gimbal frame oscilla'-^ 
tory rotation caused by beam deflection during vibration^ hardened steel 
journals were fitted into both sides of the gimbal frame at its axis of 
rotation. The journals were first put through the holes in the gimbal 
supports and then seated in the gimbal frame holes with 5/l6" shoulder 
bolts . Operation of the apparatus later indicated that the journals came 
loose in the gimbal frame holes ^ causing local heating and eventual seiz- 
ure of the frame within the gimbal supports. To correct this situation^ 
each journal was keyed to the gimbal frame ^ effectively preventing any 
relative motion between the journal and gimbal frame. In addition^ a 
journal oil supply system was Installed; oil cups were mounted on top of 
each gimbal support post with oil wicks leading down to the journal through 
holes in these support posts. To distribute the oil around the journal^ a 
circumferential l/l6" x l/l6" groove was milled in the center of the 9/16" 
diameter hole of each gimbal support post. A 10 per cent colloidal dis- 
persion of molybdenum-disulfide in oil was used as the lubricant. 

Each gimbal assembly was secured to an alignment platform by four 
3/8” bolts. Since the nodal position was located by calculation (see 
Appendix A) ^ small errors in gimbal frame moimtlng position did exist. 

Thus^ some vertical force was applied to the gimbal supports. A rigid 
connection of the gimbal assembly to the alignment platform would have 
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introduced a significant amount of damping. Therefore^ provision for 
small relative motion between the gimbal assembly and the alignment plat- 
form was made by placing a 1/32” sheet of neoprene rubber between them. 

The alignment platform was designed to insure accurate alignment of 
the resonant beam with respect to the vibration exciter driver axis. Ac- 
curate positioning of the resonant beam was necessary so that no lateral 
force component could be applied to the exciter driver coil through the 
exciter extension. Clearance between the exciter driver and field coils 
was such that any misalignment of the exciter extension could cause the 
driver coil to come in damaging contact with the field coil. The align- 
ment platform^ machined from a standard 8” structural channel section^ 
was slotted along the top center-line to receive the pad projections of 
the gimbal supports. This arrangement prevented any lateral motion of the 
resonant beam. To facilitate moving the gimbal assemblies to the resonant 
beam node positions four slots ^ 16” in lengthy were machined in the align- 
ment platform. 

Small stainless steel bench marks were located on the alignment 
platform on either side of the central 1 3/8" diameter hole for the alumi- 
num exciter extension. These bench marks served as a convenient aid in 
locating the gimbal assemblies at equal distances from the beam center. 

A special brass alignment plug was made to fit into the 1 3/8" 
diameter hole through which the exciter extension was centered. The 
alignment plug had its center drilled out to the exact diameter of the 
exciter extension. Thus^ the installation position of the alignment 
platform on the vibration exciter was accurately determined when the 
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alignment plug, with the exciter extension through its center, would just 
slide into the 1 3 / 8 " diameter hole in the center of the alignment platform. 

The alignment platform was attached to the vibration exciter by 
quadrants and brackets in an arrangement which permitted rotation of the 
entire exciter-alignment platform portion of the test stand relative to 
the exciter base. In this manner the direction of vibration could be 
varied 0 to 90° relative to the direction of gravity. The brackets were 
bolted to the top face of the exciter casting, and the quadrants were 
fitted to the same trunnions which normally connect the exciter to its 
base . 

The exciter base was in turn bolted to a 2' section of l4" x l4 l/2" 
wide-flange beam which served as the base support for the entire vibration 
test stand, as shown in the photograph. Fig. 10. Adjustable legs were 
installed at each corner of the wide-flange base to provide a convenient 
means of leveling the test stand. These legs supported the base through 
i*-" square pads of heavy rubber which prevented vibration transmission to 
the test room floor. 

The yoke which held the test cylinder at mid-span of the resonant 
beam was a U-shaped, finned structure which combined high rigidity with 
minimum weight. The base plate of the test cylinder yoke was held to the 
resonant beam by means of a steel clamping plate. The resonant beam was 
firmly clamped between these two plates held together by four \/h" steel 
bolts. A calculation of stresses showed that the maximum stresses in these 
bolts during vibration were at all times well below the allowable limit. 

To prevent any corrosive or cutting action between the plates and the 



i 







I 



I 




22 . 



resonant beam^ thin pieces of Teflon sheet were placed between the plates 
and the beam. 

The lead wires of the test cylinder were fastened to the test cyl- 
inder in the manner described in section 2-b of this chapter. These 
wires were brought down along the test cylinder yoke and taped to the 
resonant beam as seen in Fig. 11. The lead wires were held firmly along 
the yoke by taping them to a stiff U-shaped wire which was anchored to 
the yoke at both ends and at the center of the resonant beam as can be 
seen in the close-up' photography Fig. 11. The wires were then taped 
along the side of the beam and led off the beam near the nodes to cable 
clamps fitted with rubber grommets on the giiribal support posts. This 
arrangement prevented any large relative motion between adjacent parts 
of the wires y and lead wire breakage was avoided completely. 

1-d. Vibration exciter and exciter power supply 

The external driving force for vibration at resonance was supplied 
by a MB Electronics IVfodel GIF electro dynamic vibration exciter. The 
sinusoidal force output of this exciter was 50 maximum within a freq- 
uency range of 6 to 350 cps. Driver coil motion was limited to a double 
amplitude of l/2”; an over-travel limit switch was provided to open the 
driver coil circuit in case this limit was exceeded. The force generated 
was directly proportional to the value of alternating current in the 
driver coil for constant direct current input to the field coil. The 
force rating of the exciter with 1.75 amperes input to the field coil 
was 10 lbs vector force for each 1.0 ampere r.m.s. of alternating current 
applied to the driver coil. The maximum driver coil current was 5*0 
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The power supply for the vibration exciter was the associated MB 
Electronics Model TII203I Power Supply. This power console contained an 
audio -oscillator with amplifier which supplied a sinusoidal wave form 
current to the driver coil and included built-in impedance matching con- 
trols . The field coil supply in the console was a rectifier which pro- 
vided direct current up to the rated 1.75 amperes for the exciter field 
coil. 



l-»e. Amplitude and frequency measuring apparatus 

The vibration amplitude of the test cylinder was measured optically 
by sighting through a cathetometer on a specially made target. The tar- 
get consisted of a short piece of 0.001” diameter wire mounted on a 
black background. When the target was illuminated by an intense spotlight^ 
the image visible through the cathetometer microscope was a thin bright 
line corresponding to the sharply illuminated crest of the wire. When 
the cylinder was put into vibration^ this image became a bright distinct 
band with a width equal to the double amplitude of vibration. The ac- 
curacy of amplitude measurement was better than 0.002". 

Frequency measurement was obtained by feeding a signal from the 
audio -os dilator in the exciter driver coil circuit into a Hewlett Packard 
Model 52IA Electronic Counter. The counter frequency reading was always 
in good agreement with the calibrated dial of the audio-oscillator; the 
accuracy of the frequency measurement was within ^ 1 per cent. 
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^ • Heat Transfer Apparatus 

2 -a. General 

The heat transfer apparatus consisted of a test cylinder, a power 
supply to heat the test cylinder and a system of thermocouples installed 
in the test cylinder for surface temperature measurement. Special pro- 
blems were involved in the design and construction of this apparatus due 
to the high inertia loads and temperatures imposed on the test cylinder 
during long periods of vibration. 

Test cylinders previously used in other experiments involving heat 
transfer from vibrating cylinders have employed guard heaters on either 
side of the central test section heater in order to prevent axial heat 
conduction away from the test section. The complexity of this arrange- 
ment — that is, the use of three interconnecting sections, each with its 
own heating element and thermocouples — has resulted in premature fail- 
ure by fatigue . 

In order to avoid fatigue failure, it was decided to reduce the 
complexity of construction of the test cylinder. This was accomplished 
by placing only one heater within the test cylinder. Since guard heaters 
were not used to prevent axial heat conduction away from the test sec- 
tion, a metering device was installed at each end of the cylinder to ac- 
curately measure the axial conduction heat losses to the test cylinder 
supports . 

The test cylinder thermocouple system was designed to permit re- 
placement of the individual thermocouples should a failure occur. A 
detailed description of the heat cylinder and its components is given 



below. 









•'1 




J 



25 . 



2 -b . Test cy linde r 

The test cylinder,, shown in detail in Fig. 5^ consisted of a central 
electrically heated test section which extended between end supports that 
were designed to fit in the test cylinder yoke clamps. The test cylinder^ 
7 / 8 " in diameter^ had an overall length of 8 l/ 2 " and weighed I .06 lbs. 

The central section^ called the test section^ was 5*20" long. 

The electric heater^ installed within the test section^ was made by 
winding gage No. 32 Nichrome V wire (nominal resistance,, 10.26 ohms per 
foot) within a helical groove machined in the surface of a I/ 2 ” diameter 
ceramic base. The helical groove had a pitch of 20 turns per inch and 
served to insulate the heater wire. The heater leads consisted of two 
pieces of O. 085 " diameter drill rod which were centered at the ends of 
the ceramic base . The heater wires were attached to these drill rod leads 
with a special medium temperature solder. Regular half-and-half solder 
could not be used because of the high temperatures involved; silver solder 
was impractical because brazing temperatures in excess of lOOO^^F would 
bum the No. 32 Nichrome V heater wire. The solder used was "TEC" solder^ 
manufactured by the Handy & Hannan Brazing Products Company. This solder, 
composed of 5 cent silver and 95 cent cadmium, has a melting 
point of 640^F — about 200^F above that of the available lead-tin solders 

The heater tube was closely fitted inside a copper sleeve, O.D. 5 / 8 ” 
and held in its axial position by Teflon plugs through which the drill 
rod heater leads extended as is clearly shown in Fig. 5* This assembly, 
with an O.D. of 5 / 8 ", was then fitted into an aluminum sleeve with an 
O.D. of 7 / 8 " which was the outside element of the test section. The 
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copper sleeve between the ceramic heater base and the alumimim walls of 
the test section provided for a more uniform axial temperature distribu- 
tion along the test section. 

Section X-X of Fig. 5 shows the detail of the test cylinder end 
supports and the copper-Teflon-copper concentric rings which formed the 
heat meter for measuring heat loss by conduction to the test cylinder 
yoke supports. These concentric rings served as a heat meter in the 
following manner. The inner copper ring was at the temperature of the 
test cylinder end; similarly^ the outer copper ring was at the temperature 
of the test cylinder yoke support. The temperature distribution along the 
copper rings ^ a length of 1 3/8”^ “^as assumed to be constant because of 
the excellent thermal conductivity of copper. The Teflon insulating ring 
produced a temperature drop between the two copper rings . This tempera- 
ture differential was measured by a differential thermocouple made by 
soldering a gage No. 24 constantan wire between the two copper rings as 
shown in Section X-X^ Fig. The heat loss from the test cylinder by 
conduction to the yoke supports necessarily passed through these heat 
meters. This heat loss^ proportional to the temperature differential^ 
^thm^ was measured in the manner described in Chapter V. 

The copper and Teflon rings were assembled on the test cylinder ends 
by press-fitting them firmly in place. Each heat meter thus formed a sup- 
port at the test cylinder end which was held in place by the test cylinder 
yoke clamp. Fig. 11 provides a clear view of the final assembly. 

The ends of the test cylinder were sealed by Teflon terminal plates 
which were fitted into each end of the alumimmi sleeve a distance of 
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1 l/8". Small Teflon washers, l/l6" thick, separated the Teflon terminal 
plates and the Teflon plugs in the ends of the copper sleeve surrounding 
the heater tube. The outer 2 l/l6" diameter flange of each terminal 
plate, shown in Fig. 11, was the surface on which all electrical lead 
connections were made. The heater lead drill rods were supported by l/2" 
long stainless steel plugs which were threaded into each Teflon terminal 
plate. The heater lead connections from the power supply were secured to 
these plugs by means of special nuts through which set screws were placed 
so as to bear against the drill rod heater leads, thus insuring good 
electrical contact. 

The measurement of test cylinder surface temperature was accomplished 
by means of thermocouples placed inside the aluminum sleeve walls. The 
test section was provided with six thermocouples arranged as indicated in 
Fig. 5 Sections B-B through F-F. Two additional thermocouples were 
installed under the inner copper ring of each sandwich at Sections A-A 
and G-G. Thus, these two thermocouples measured the temperature of the 
inside copper ring within a close approximation. 

The thermocouples were made of gage No. 30 dual conductor copper 
constantan wire encased in fiberglass insulation. Each thermocouple junc- 
tion was made by removing the insulation from both wires along about 1" 
at the mid-span of a piece of thermocouple wire approximately 12" long. 

The two wires at the clean section were then twisted tightly together. 

Next, one end of the insulated portion of the wire was carefully pulled 
through a section of O.O65" diameter stainless steel tubing which had been 
copper-plated on its outside surface to a diameter of O.068". When the 



A 



» 



I 

f 






( 



28. 



twisted Junction was pulled just inside this tuhe^ Handy & Hannan TEC 
solder was used to encapsulate the junction at the tip of the tube. Each 
tube was then carefully fitted into its individual hole in the aluminum 
sleeve wall. These holes^ O.O 68 " in diameter^ received the thermocouple 
tubes in a tight fit; the soft copper-plating on the surface of the tubes 
provided excellent thermal contact between the tubes and the test cylin- 
der wall. 

The lead wires of the thermocouples and the two copper wires used 
for each heat meter described above were closely fitted through holes 
drilled in each Teflon terminal plate. From these holes the wires were 
led to individual terminal affixed to the terminal plate by No. 3 stain- 
less steel screws which were threaded into the Teflon. Each wire end 
was formed into a spiral similar in form to the hair-spring of a watch. 

The spirals were secured with small copper strips held by the No. 3 
screws. The spiral end prevented any wire from crossing itself and break- 
ing from the contact when under intense vibration. The lead wires lead- 
ing down the test cylinder yoke (described in Section 1-c of this chapter) 
were attached in a similar manner. Thus^ each wire terminal consisted 
of spiral end from the thermocouple junction, a copper strip, the spiral 
end from the lead wire, a second copper strip and, finally, two No. 3 
screws holding the terminal to the Teflon terminal plate. To prevent 
these screws from coming loose during vibration, their slots "^ere aligned 
and small pieces of wire were soldered across the slots, preventing ro- 
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The extreme care used in the installation and wiring of the thermo- 
couples proved worthwhile; not one thermocouple failure occurred during 
the test program. However, an attractive feature of this thermocouple 
installation method was that an individual thermocouple could be removed 
and replaced in the event of failure . 

The surface of the test section was polished to a mirror finish to 
insure a low thermal emissivity. A low value of emissivity made the 
radiation heat loss from the test section small compared to the heat 
transferred by convection. Consequently, errors made in the calculation 
of radiation losses did not introduce appreciable errors in the calcula- 
tion of the convective heat transfer coefficient. The emissivity, as 
determined by comparison with a similar heat transfer test cylinder of 
known emissivity, was e = 0.0^5 . 

2-c. Power supply 

The power supply to the test cylinder heater was drawn from a ^00 
watt constant voltage regulator. Two variacs, placed in series with the 
voltage regulator output, served for the coarse adjustment of the voltage 
across the heater. A fine adjustment was provided by a one ohm variable 
resistor in series with the heater. The variacs and the variable resis- 
tor were mounted on a panel together with an ammeter and voltmeter, and 
appropriate switches, fuses and indicator lights. 

The power supplied to the test cylinder heater was measured con- 
tinuously by an external Weston Model ^32 double -current range wattmeter. 
The meter was connected in such a manner that the meter indication 
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included only its potential circuit losses. On low scale, used to measure 
power below 30 watts, the internal resistance of the potential circuit was 
4083 ohms, as measured by a high precision portable Wheatstone bridge. On 
high scale, used to measure power above 30 watts, this resistance was 8180 
ohms. The resistance of the heater was 9O.O7 ohms; the resistance of the 
heater leads was 0.17 ohms. Thus, with the meter on low scale, the po- 
tential and lead circuit loss was 2.l6 per cent of the wattmeter Indication 
on high scale the loss was I.09 per cent of the wattmeter indication. 

2-d. Temperature measuring apparatus 

The eight thermocouples mounted in the test cylinder walls and the 
two additional thermocouples used to measure ambient temperature in the 
test room were connected through a Minneapolis -Honeywell thermocouple 
selector switch to a common cold junction and a portable Rubicon precision 
potentiometer. The cold jmction was placed in a thermos bottle filled 
with water and ice. 

The differential thermocouples of the two heat meters at the test 
cylinder ends were connected to the Rubicon precision potentiometer 
through a double-pole, double-throw knife switch. The thermocouple 
selector switch was put on open circuit when it was desired to obtain 
the heat meter temperature differential at each end of the test cylinder. 

Standard tables were used to convert the potentiometer readings to 
temperatures. It was considered that errors associated with determining 
absolute temperatures and temperature differences with this system were 
within 0.1°F. 
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CHAPTER IV 

EXPERIMENTAL PROCEDURE 

The experimental program was carried out in two phases. The first 
phase, called the calibration phase, was a series of tests performed with 
no vibration in order to determine the heat loss by conduction, Q^, to 
the test cylinder supports. The second phase, called the vibration phase, 
involved the collection of heat transfer data with the test cylinder in 
horizontal vibration. 

The calibration phase was carried out with the test stand and test 
cylinder in the horizontal position as shown in the photograph. Fig. 10. 
The test section heat input was set at integral multiples of 5 vatts 
from 5 to 50 watts. Two additional tests were conducted at wattmeter 
settings of 60 and 70 watts . The rated accuracy of the wattmeter was 
0.5 per cent of full scale deflection; hence, the maximiom possible error 
of the wattmeter reading occurred at the 5 "w-att level and was 3*5 por 
cent. The possible error at the 70 watt level was 0.5^ per cent. 

Since the heat input to the test cylinder was kept constant during 
a given test, the steady state temperatures of the test cylinder were 
approached exponentially in time. The criterion established for terminat- 
ing the calibration phase tests was a maximum change in any individual 
test cylinder thermocouple reading of 0.2^F and a maximim change in heat 
meter temperature differential, Atj^, of O.l^F over a period of one- 
half hour. In practice, a steady state temperature was generally achieved 
with less variation than the above requirements. The average calibration 
test run extended for a period of approximately four hours . 
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The ambient temperature, t . was given by a thermocouple placed 
about 15” below the test section. A second thermocouple^ suspended in 
the test room about 3 above and behind the test cylinder was monitored^ 
and generally read within 0.2^F of the thermocouple below the test sec- 
tion. 

The surface temperature of the test section was defined as the 
average of the six temperature readings from the six thermocouples in the 
test section locations indicated in Chapter III, Section 2-be As expect- 
ed, the axial distribution of the individual test section teirperature 
readings showed a variation dependent on heater power level. The variation 
between the center and end individual test section temperatiires ranged 
from 4.0^F at the 5 vatt level to 40c3°F at the 70 watt level. The tem- 
perature potential. At, based on average surface temperature, t^^ from 
the six test section thermocouple readings was 25»6^F for the 5 watt level 
and 273 -9°F for the 70 watt level. Thus, considering the worst condition - 
that the average tg was in error by one-half of the axial temperature 
range (2^F for the 5 watt case and 20.15^F for the 70 watt case) - the 
percentage error in At would be 7*8l per cent at the 5 watt level and 
7.38 per cent at the 70 watt level. These errors, of course, represent 
the unlikely extreme; in actual test, the measured average test cylinder 
surface temperature was estimated to be within 2 per cent of the true 
mean surface temperature . 

The temperature differential across the heat meters at the test 
cylinder supports was provided by the differential thermocouples (described 
in Chapter III, Section 2-^b). The temperature level from which the tem- 
perature differential was measured was provided by the thermocouple 
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immediately beneath each heat meter in the test cylinder surface (at Sec- 
tions A-A and G-G^ Fig* 5 )<» The left and right heat meter temperature dif- 
ferentials were always within 2^F of each other. This variation was due 
to a small difference in heat flow across the left and right test cylinder 
heat meters. In practice^ the average of the two heat meter temperature 
differentials was taken as the test cylinder end temperature difference, 

At, .As stated in Chapter III, Section 2-d, it was considered that 
hm 

the errors associated with determining absolute temperatures and tempera- 
ture differences with this system were within ± O.l^F. 

The results of twelve calibration phase test runs are given in 
Table 1 and plotted in Fig. I 3 . The method used in calculating these 
results is presented in Chapter VI. 

The vibration phase of the experimental program was oriented toward 
the major objective of the entire investigation; namely, the investigation 
of the influence of horizontal mechanically induced simple harmonic vibra- 
tions upon the heat-transfer rate by free convection from a heated hori- 
zontal cylinder to air. Specifically, the tests were so arranged as to 
provide data which could be compared directly with the sound method data 
obtained by Fand and Kaye as explained in Chapter I. 

The three principal experimental variables were as follows: 

a. temperature potential, ^t 

b. intensity of vibration, af 

c. heat input to the heat transfer cylinder test section, 

In order to provide a direct comparison with the sound method experiments 
of Fand and Kaye in a small number of test runs, it was decided to restrict 
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the independent variable At to one value^ namely^ At = 150 ^Fo Since^ 
in practice^ it would be extremely difficult to obtain an exact steady 
state value of = 150 *^? for each test^ the allowable variation from 
this temperature potential was established as ± It was also nec- 

essary to vary af within a range of about 0.2 to 1,1 ft/sec to cover the 
same range of vibration intensity as was done in the sound method experi- 
ments . 

The vibration tests we\re all made at the same frequency^ f = lOk cps^ 
the resonant frequency of the mechanical system. At the start of each 
test^ the vibration exciter was tuned to this resonant frequency. Tuning 
was achieved by applying the principle that for a coupled electro- 
mechanical system with a constant electrical energy inputs the current 
becomes a minimum at resonance . An additional aid in determining reson- 
ance was the fact that the maximum apparatus noise levels while not large ^ 
coincided with the resonant condition. 

After obtaining vibration at resonance^ the desired value of af was 
acquired by adjusting the output level of either the audio -oscillator or 
its amplifier or both. Also^ at lower vibration intensities it was found 
convenient to limit the value of the direct current output to the exciter 
field coil to an arbitrary value less than the 1.75 ampere rating. This 
procedure^ in effect^ added sensitivity to the gain controls of the 
audio -os dilator and amplifier and thereby aided in establishing the de- 
sired amplitude. These gain controls were adjusted in trial and error 
fashion until the double amplitude of vibration was at the value required 
to give the proper af value for the test run. This procedure^ initially 
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tedious, was usually accomplished within 10 minutes after some experience 
had been acquired » The vibration amplitude was checked at 5 minute inter- 
vals; minor correction was occasionally necessary. 

The requirement of producing a steady state value of At within l^F 
of 150^F was more difficult to achieve. Because the steady state heat 
input to the test cylinder with vibration was unknown, its value had to 
be determined in each test run. The power supplied to the test section 
which resulted in a 150 ^F temperature potential in the absence of vibra- 
tions was determined to be about 3^ watts from the results of the calibra- 
tion tests. It was assumed that the power level needed to produce the 
same temperature potential with the test cylinder in vibration was greater. 
The test cylinder was heated for about two hours prior to vibration at a 
power level somewhat greater than 3^ watts . After the desired vibration 
intensity was established, a complete set of temperature readings was 
taken at 10 minute intervals. The temperature potential existing at the 
time of temperature measurement was then quickly computed, and the power 
level was adjusted to bring At closer to 150^F. This process was re- 
peated until no change in A t greater than O.l^F in a 10 minute period 
was detected* At this time the watt meter reading was recorded as the 
steady state value, and the cathetometer readings of vibration amplitude 
were recorded. 

The dToration of the individual vibration tests varied between two 
and three hours. The results of the nine tests with vibration are re- 
corded in Table 2 and plotted in Fig. l4. The method used in calculating 
these results is presented in Chapter V. 



I 




9 ^ 



i 















36 . 



CHAPTER V 

CALCULATION PROCEDURE AND RESULTS 

This chapter contains a description of the methods and basic caleula* 
tion procedure used to evaluate^ and interpret the experimental data^ and 
a description of the results obtained. 

The calibration phase tests of the experimental procedure provided 
data with which the axial heat losses from the test cylinder could be de- 
termined as a function of average heat meter temperature differential^ ^^hm” 
During the tests in the absence of vibration, the heat input to the test 
section, was dissipated to the system environment by the following modes 
of heat transmission: 

a. Heat loss from the test section by radiation, 

b. Heat loss from the test section by convection (in the absence 
of vibration), 

c. Heat loss to the test cylinder supports, 

d. Heat loss by convection from the end terminal plates of the 

test cylinder and heat loss by conduction through the test 
cylinder lead wires (in the absence of vibration)^ called 

the unaccoxmtable heat loss. 

The test section heat losses by radiation and convection, and 
Qco^ were subject to accura;te calculation as Will be described belbw. 
Calculation of the heat loss to the test cylin4er supports, Q^, could not 
be accurately calculated, however, because: (l) the thermal conductivity 

of the heat meter Teflon ring was not known within the required accuracy. 
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and ( 2 ) the contact coefficient of the copper-Teflon heat meter junctions 
a function of test cylinder support clamping pressure^ surface roughness 
and the presence of oxide films — could not be determined. Further^ cal° 
culations for using the best available value of thermal conductivity 
for Teflon^ gave results that did not agree well with the empirical deter- 
mijiation of as described below. Similarly^ could not be calculated 

because of the complicated geometry of the outside surface of the Teflon 
plates . 

The sum^(0^ + could be determined^ however^ through application 

of the following heat flow continuity equation: 



In the calibration phase tests ^ (Q^ + ^uo) determined as a function 
of At^, The results of these tests are plotted in Fig. 13« The origi- 
nal experimental data for the calibration tests and the associated cal- 
culated quantities are listed in Table 1. The experimentally determined 
curve of (Q^ + Q^q) ^^^sus iA'^hm very nearly a straight line which 
indicates that the predominate heat loss at the test cylinder ends was by 
conduction^ i.e.^ + Qqq* This was' the expected result since the out- 

side surface temperature of the Teflon terminal plates was much less than 
the temperature of the test cylinder. Also^ the heat loss by conduction 
through the lead wires was negligible compared to Q^. 

The heat transferred by convection in the presence of vibrations^ 
was determined from the heat continuity equation: 






( 6 ) 



*^CV Qv " ^ ^v) 



(T) 
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where is the unaccountable heat loss in the presence of vibration. It 
is recognized that was somewhat greater than Q^q increase 

in convective heat transfer from the terminal plates in the presence of 
vibrations; however^ this difference, (^uv " ^o) assumed negligible 
since the uinaccountable heat loss itself was small. Accordingly, in the 
tests with vibration, was computed by: 

^cv ” ^uo) 



where determined from Fig. 13 . 

The heat transferred by convection from the test section in the ab- 
sence of vibration, ’W'as determined in the following manner. The de- 

fining equation of the heat-transfer coefficient in the absence of vibra- 
tion, iiQ, is: 

Q 

CO 2 

= 3.412 : — , Btu/hr-ft -°F (^l) 

° A At 

where 

2 

A = test section heat transfer area, 0.0993 
= heat transferred by convection, watts 
At = temperature potential, ^F 
3.412 is the conversion factor from watts to Btu/hr 



(9) 

The experimentally determined equation for h^ of a cylinder of 



diameter is: 



hQ = 0.225 (■^^) ^ Btu/hr-ft^ -°F 

o 



( 10 ) 
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Equation (lO) was the result of a heat-transfer study with a cylinder of 
identical diameter il/Q”), length and end support method and surface thermal 
emissivity as the test cylinder used in this project. 

Equating Equations (9) and (lO) with = 0.0729 ft and solving for 



gives: 



A 

= 0.0143 A t , watts 



( 11 ) 



The maxim'um experimental error in the determination of h^ by Equa- 
tion (lO) was reported to be 2.6 per cent.^^^ The maximum possible experi- 
mental error in the calculated value of was determined by logarithmic 
differentiation of Equation (9): 

dh dQ dA d(^t) 

O CO / . 

»o ' «co ' A ' 



Equation (l2) was written in finite difference form; hence, the maximum 
possible percentage error in was the sum of the maximum possible per- 
centage errors of the quantities involved in its calculation. The possible 
error in h^ was 2.6 per cent; the possible error in the calculation of 
surface area. A, was less than 0.2 per cent. The maximixm error in ^t was 
estimated to be 2 per cent (see Chapter IV). Thus, the maximum possible 
error in the value of was 4.8 per cent. 

The heat transferred by radiation from the test section, was 

independent of the presence of vibrations. Q^, as a function of At, was 
determined from the graph. Fig. 12, based on the equation: 



Qr = 




0.2931 CTe A 



, watts ( 13 ) 



imp 

pm. 








'■ 1 ' 



I 



ho. 



where 



o~ 

e 

A 

A ^ 
0.2931 



Stefan-Boltzman constant, 0.1713x10”® Btu/ft^-hr-°R^ 

surface thermal emlssivity of the test section, 0.045 

2 

test section heat transfer area, 0.0993 ft 
mean ambient temperature, 5^0°R 
temperature potential, °F 
conversion factor from Btu/hr to watts 



The errors associated with determining Q^, which were introduced by the 
fact that the actual ambient temperature was in most cases slightly dif- 
ferent from the mean ambient temperature, were less than 4 per cent. 

These errors were negligible because itself was small compared to the 
total heat transferred from the test section. 

The heat transfer coefficient in the presence of vibration, h^, was 
calculated for each test of the vibration phase by using the heat trans- 
ferred by convection, for the test as determined by Equation (8). 

Thus : 

^cv 

\ = 3.412 f watts (l4) 

The maximimi possible experimental error in the calculated value of 
h^ was determined by writing the logarithmic differential equation of (l4) 
in finite difference form as in the case for above. The maximum pos- 
sible error involved in obtaining was 3-0 per cent, as obtained by 
adding the maximum errors possible in and + Q^o)* Errors in 
were less than 0.1 per cent and were neglected. Again the errors in sur- 
face area. A, and ^ t were taken as 0.2 and 2 per cent, respectively. 
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Thus^ the maximum possible experimental error in the determined value of 
h^ was 5*2 per cent. 

The experimental data and the calculated results of nine tests with 
vibration are recorded in Table 2. The range of vibration intensity em- 
ployed was from 0.213 to 1.109 ft/sec; the temperature potential in all 
these tests was within ± l^F of 150^F. The value of h^ varied from 
1.803 Btu/hr-ft^-^F at the lowest af to 4.323 Btu/hr-ft^-^F at the highest 
af . The value of hQ at ^ t = 150^F was computed to be 1.71? Btu/ 

hr-ft^-^F by Equation (lO). Thus^ percentage increases in the heat-transfer 
coefficient from 5«0 to I5I.8 per cent were achieved. 

In order to provide a direct comparison of these results with the 
sound method investigation of Fand and Kaye^^^^ (discussed in Chapter l)^ 
the variation of h^ with af at A t = 150^F for a stationary horizontal 

cylinder in the presence of a high intensity sound field of frequency 
f = 1496 cps was obtained from Fig. 6l of Reference 10 by cross -plotting. 

The data points of the vibration tests of the present investigation were 
then plotted on the same graph^, Fig. l4. The maximum deviation of any 
single data point from the sound method curve in Fig. l4 is 4.6 per cent 
as compared to the maximum possible experimental error of 5*2 per cent in 
obtaining the data points. Fand and Kaye arrived at a correlation equa- 
tion which expressed the dependence of h^ on At and af with a high degree 
of consistency above an af of about O.7 for sound frequencies less than 
2000 cps. For convenience this correlation equation is restated here: 



h 

V 



^ t (af 



1/3 



0.722 



( 1 ) 



1 1 



t 
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The maxiffluni deviation of h^ obtained in this investigation for the three 
values of af greater than O .7 from the computed value h^ given by Equa- 
tion ( 1 ) was 7-3 cent. The average deviation was 2.1 per cent. 

A comparison of the results of this investigation with the results 

of the study of heat-transfer rate variation with vertical mechanical 

( 8 ) , 

vibrations (discussed in Chapter I) was also made. A curve of h^ ver- 
sus af for At = 150^F was obtained from Fig. 12 of Reference 8 and is 
plotted in Fig. I 5 of the present report. The data points of the vibra- 
tion phase of this investigation were plotted on Fig. I 5 to facilitate com- 
parison. Inspection of Fig. I 5 revealed that the trend of h^ increase 
relative to af was quite similar for the two investigations. In particu- 
lar^ at values of af above about O .3 ft/sec ^ the value of h^ is somewhat 
less for the horizontal vibration case^ while the slopes correspond closely. 
This result will be discussed further in Chapter VI. 

It was considered desirable to obtain some visual evidence of the 
type of boundary-layer flow that existed around the heated test cylinder 
during vibration. Accordingly^ a crude flow- visualization examination was 
performed^ using a lighted cigarette held beneath the cylinder as a smoke 
source. The results of this visual observation clearly indicated that 
laminar flow persisted in the region near the lower stagnation point 
throughout a range of af from 0 to about 1.1 ft/sec. Also^ some evidence 
of vortex formation above the cylinder was noted; this vortex flow is 
characteristic of thermo acoustic streaming as observed by Fand and Kaye 
in their sound studies. Although the results of the elementary visualiza- 
tion technique employed here require further substantiation^ the study 
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did indicate conclusively that the fluid flow mechanism in the region of 



the lower stagnation point of the cylinder was not vih rationally induced 
turbulence as reported for the vertical vibrations case^^^ 
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CHAPTER VI 

DISCUSSION OF RESULTS 

The results of the present experimental investigation reaffirm the 
conclusions of previous investigations that the influence of vibrations 
on heat transfer is significant. The plots of h^ versus af (Figs. l4 and 
15 ) clearly indicate the general effect of the Intensity of vibration on 
the heat-transfer rate from a heated horizontal cylinder in horizontal 
transverse vibration in air. It was noted that the change in h^ at vibra- 
tion Intensities less than af * O .3 ft/sec is not appreciable. Above this 
value of vibration intensity, the heat-transfer coefficient increases 
rapidly with increasing af. The value of af = O .3 ft/sec, where vibra- 
tions take on major significance with regard to the heat-transfer rate 
is called the critical intensity. The critical intensity, af = O .3 ft/sec, 
observed in this investigation corresponds very nearly to the critical 
intensities reported for the sound method (Fig. l4) and also in the ver- 
tical mechanical vibrations method (Fig. 15)* 

The close agreement of the quantitative heat transfer results and 
the flow visualization observations of this investigation^ with those of 
the soimd method for comparable vibration vector orientation provides 
evidence which indicates that the physical mechanism of the heat transfer 
for the two cases is identical. However, the quantitative results, when 
compared to the vertical vibrations method (Fig. 15)^ show that the 
slopes of hy. versiis af are very nearly the same for these two cases. The 
fact that the trend of the data obtained in this investigation is similar 
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to the results of the vertical vibrations method (Fig. I 5 ) is not con- 
sidered evidence that the two modes of vibration are equivalent with 
respect to heat transfer; for, as has been explained in Chapter V, the 
mechanism of interaction between vibrations and free convection in the 
case of vertical vibration is vibrationally induced turbulence, whereas 
there is evidence that the mechanism of interaction for both the sound 
method and the present investigation is thermoacoustic streaming. 

The slope of h versus af for af greater than about 0.9 for the 

V 

sound method is 2/3 as Equation (l) indicates. The corresponding slope 

(9) 

for the vertical vibrations method has been reported to be 1 . Addi- 

tional heat transfer data is required for values of af greater than O .9 
to determine conclusively the value of this slope for the horizontal 
vibrations method. 

The ratio of amplitude of vibration to cylinder diameter, a/D^, 
is a controlling variable in all fluid-mechanical phenomena involving 
vibrations. The a/D^ ratio for this investigation at a typical value 
of af = 0.9 ft/sec and f = 104 cps was 0 . 119 ; the a/D^ ratio for the 
sound method of Fand and Kaye with af = O .9 ft/sec at f = 1496 cps was 
9.63 X 10'3. Thus, the a/D^ ratios differed by a factor of 12. 3 . In 
view of the close quantitative agreement of the results of the two 
methods of investigation (Fig. l4), it appears that the influence of 
u/Dq does not have a predominant effect, at least for the range of a/D^ 
ratios indicated here . 
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CHAPTER VII 
CONCLUSIONS 



The experimental portion of this investigation of the influence of 
horizontal mechanically induced vibrations on the heat-transfer rate from 
a horizontal cylinder to air consisted of a total of 21 heat transfer 
test runs. Twelve of these test runs provided calibration data on axial 
heat loss from the heat transfer test cylinder; nine test runs were con- 
ducted with the cylinder in horizontal vibration. The ranges of the 
principal experimental variables were: 



temperature potential^ At: 
frequency of vibration^ f: 
amplitude of vibration^ a: 
intensity of vibration^ af: 



150 °F + 1 °F 
10^ cps 

0.025 to 0.128” 

0.21 to 1.11 ft/sec 



The results indicate that the influence of horizontal vibrations 
upon the heat-transfer rate by free convection to air becomes appreciable 
when the intensity of vibration^ defined as the product af^ exceeds the 
critical value of O.3 ft/sec. This critical value is approximately equal 
to the critical sound vibration intensity determined in Reference 10 and 
the critical value of af for vertical mechanical vibrations determined in 
Reference 9 * The flow near the horizontal cylinder in the presence of 
intense mechanically induced transverse horizontal vibrations was ob- 
served to be laminar near the lower stagnation point using smoke as the 
indicating medium; this type of flow is one characteristic of the flow 
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near the lover stagnation point of a cylinder in the presence of intense 
acoustically induced transverse horizontal vibrations. The observed 
laminar flow near the lover stagnation point vas quite different from the 
vibrationally induced turbulent flov characteristic of horizontal cylin- 
ders in transverse vertical vibration. 

Based on the quantitative heat transfer results and the elementary 
flov study of this investigation, it appears that the physical mechanism 
of heat transfer for the sound and mechanical vibrations methods, both 
vith geometry as indicated in Fig. 1-a, is identical. Also, the com- 
parison of these tvo methods indicates that the amplitude of vibration 
to cylinder diameter ratio, a/D^, does not have a predominant effect for 
the range of a/D^ ratios considered. 
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CHAPTER VIII 
RECOMMENDATIONS 



This chapter presents recommendations for future investigations in 
the field of heat transfer from vibrating cylinders based on the results 
and experience gained from the present project. Six specific suggestions 
are made; the first five of these relate to experimental work that can be 
conducted with the apparatus designed and constructed for this project. 
Indeed, this apparatus was designed to include features whereby adaptation 
to other investigations will require minor or no modification. The list 
of suggested investigations with a few relevant remarks is given below. 

(1) Extension of the range of vibration intensity to values of af 
greater than 1.1 ft/sec. This is necessary, as indicated in Chapter VI, 
in order to determine accurately the slope of h^ versus af in the higher 
af range. The existing apparatus has been tested at af = I.3 ft/sec 
with about 35 exciter driving force. It is considered that an af of 
at least I .5 ft/sec can be achieved at the maximum exciter output of 50 
lbs force , 

( 2 ) Flow visualization studies. It would be of particular value to 
obtain photographic evidence of the flow pattern around the test cylinder. 

(3) Change in vibration direction relative to the direction of the 
force of gravity. Since the flow character aroimd the cylinder has been 
shown to be different for the vertical and horizontal vibration cases, it 
would be of value to investigate vibration angles between these two ex- 
tremes . This can be done with the existing apparatus as described in 
Chapter III, Section 1. 
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( 4 ) An investigation of velocity and temperature distributions 
near the cylinder. It is considered that hot-wire anemometry techniques 
could be used to determine detailed information concerning the heat 
transfer mechanism near the cylinder by an investigation of this type. 

(5) Different cylinder diameters. The quantitative determination 
of the influence of cylinder diameter is important. New test cylinders 
and some modification of the test cylinder yoke would be required. 

(6) Analysis. It is considered that^ on the basis of the data 
provided herein and from other investigations in the fields an analysis 
should be attempted which would provide a mathematical foundation for 



the phenomena observed. 
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FIG. I ORIENTATION OF VIBRATION VECTOR RELATIVE 
TO DIRECTION OF GRAVITY FOR (a ) HORIZONTAL 
TRANSVERSE VIBRATIONS AND (b) VERTICAL 
VIBRATIONS OF A HORIZONTAL CYLINDER 
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NOTES : 

1. DIMENSIONS IN BRACKETS INDICATE NOMINAL 
SIZE OF MATERIAL, AS DELIVERED. 

2. GRIND BOTTOM 5/16" EITHER SIDE OF 
SLOT along LENGTH OF SLOT. 

3. INDICATES POSITION OF QUADRANT 
SUPPORT PLATE CONTACT ON CHANNEL 
BOTTOM. GRIND THIS AREA FLAT. 

4. DRILL OOWEL PIN HOLES IN PARTS 
I a 2 AT THE SAME TIME. 

5. DEFINES LOCATION (UNDER) OF SPOT 
FACE TAPS ON EXCITER CASTING TOP 
AND THE LENGTH-WISE DIMENSION 
FROM WHICH THE 4 FLANGE TAPS 
SHOULD BE LOCATED. 

6. PART 3a TO HAVE LIGHT PRESS FIT 
INTO PART 3e a SLIDING FIT THROUGH 
PART 4. 






I ALIGNMENT PLATFORM 




C‘ 



0 I2S' 



3c. BEAM SLIPPER FOR l.43x 1.50 BEAM 

3e GIMBAL FRAME 



SCRIBE 
TOP t 



NO. 30 DRILL, CSK FOR. 
I/B" FLAT HEAD SCREVV, 




2 BENCH MARK 



U-7/eJ 



SET (2) l/B' OIA. OOWEL PINS, 
PRESSEO INTO PART 2, FOR 
SNUG FIT INTO PART I. 



DESCRIPTION OF PART 


NO. 
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MATERIAL 


REMARKS 


1. ALIGNMENT PLATFORM 


1 


STRUCTURAL STL 


BljIB.S 4-0" 


2.BENCH MARK 


2 


303 STAINLESS 




3a. GIMBAL JOURNAL 


4 


1020 ST. H.R. 




3b. BEAM CLAMP PLATE 
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303 STAINLESS 




3c. BEAM SLIPPER 


2 


303 STAINLESS 


FOR I.43XI.S0BEAM 


3d. BEAM SLIPPER 


2 


303 STAINLESS 


FOR I.OOxl.SO BEAM 


3e. GIMBAL FRAME 
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2024-T4 AL. 






VIBRATION TESl 


STAND 


DETAIL DRAWINGS OF 
PART 1, 2 a 3 


SHEET 2 OF 3 
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DRILL (2» 9/16" OIA 
HOLES TO TAKE 
PART -30 IN SLIDING FIT 



SEE NOTE 4 
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jl/2" 



4 . 000 **^^^ 
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SEE NOTE 2 
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8 "RAD. 
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KET 



NOTES 

1. PART 4. PART 3 TO FIT WITHIN POSTS WITH 
0.001" MAX. CLEARANCE 

2. PART4. PAD PROJECTION TO HAVE CLOSE, 
SLIDING FIT WITHIN SLOT OF PART I 

3. PTS 3 8 8. COMPLETE ALL WELDING 8EF0RE 
DRILLING a MACHINING 

4. GRIND LOAD SURFACE FLAT WITH THE 2 
QUADRANTS FIXED 8ACK-T0-8ACK 



DESCRIPTION OF PART 


NO • 
REO'D 


MATERIAL 


REMARKS 


4. GIM8AL SUPPORT 


2 


304 STAINLESS 


1/2" STOCK 


5. LOAD PLATE YOKE 


1 


1020 ST H. R. 


5/8" XM/2" RECT ROD 


6. LOAD PLATE 


1 


RY - ALLOY 


5/32" X 3/4" GROUND STOCK 


7. CLAMPING PLATE ASS’Y 


1 


1020 ST. H.R. 




8. QUADRANT 


2 
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1/4" PLATE 


9. BRACKET 


2 


MILD STEEL 


1/4" PLATE 


10. EXCITER EXTENSION 


1 
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VIBRATION 


TEST STAND 


DETAIL DRAWING OF PARTS 4-10 SHEET 3 OF 3 
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FIG. 5 TEST CYLINDER 
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BEAM TO CENTER LOAD MASS RATIO VS Q-e AT RESONANCE 



Fig. 6 
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NODAL DISTANCE FROM BEAM END TO HALF BEAM LENGTH RATIO VERSUS a-l 
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maximum bending stress at beam center, 
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af , f ps 

0 0.25 0.5 0.75 1.00 1.25 1.50 1.75 2.00 




MAXIMUM BENDING STRESS VS AMPLITUDE AND af FOR 
STEEL AND ALUMINUM BEAMS AT VARIOUS VALUES 
WITH A FREQUENCY OF 100 CPS 



Fig. 9 
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Fig. 10 - Vibration Test Stand, Horizontal Position 
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Fig. 11 - Close-up View of Test Cylinder 
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HEAT TRANSFERRED BY RADIATION, Q, , watts 
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HEAT TRANSFERRED BY RADIATION FROM TEST SECTION 



Fig. 12 
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HEAT TRANSFERRED FROM TEST CYLINDER ENDS, (Qd + Quo) 

VS A t hm 



Sk. 



Horizontal cylinder subjected to horizontal 

transverse acoustical vibrations, Do = 3/4", 
f = 1496 cps, At= I50°F 

O Horizontal cylinder subjected to horizontal 
transverse mechanical vibrations, Do = 7/8" 
f = 104 cps. At = I50°F 




a f , f t / sec 

FIG. 14- COMPARISON OF hv vs at AT CONSTANT At FOR 

HORIZONTAL CYLINDER IN HORIZONTAL TRANSVERSE 
MECHANICAL AND ACOUSTICAL VIBRATION 
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Horizontal cylinder subjected to vertical 

transverse vibrations, Do= 7/8", At = 150® F 

0 Horizontal cylinder subjected to horizontal 
transverse vibrations. Do =7/8", At = l50®F 



a f j f t/sec 



JO .2 .3 .5 .6 .7 .a .9 i.O 




FIG. 15 COMPARISON OF hv vs at AT CONSTANT At FOR 

HORIZONTAL CYLINDER IN HORIZONTAL AND VERTICAL 
TRANSVERSE MECHANICAL VIBRATION 
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EXPERIMENTAL AND COMPUTED VALUES OF CALIBRATION PHASE TESTS 
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APPENDIX 



I 

I 

u 



TO. 



APPENDIX A 



VIBRATION ANALYSIS 



In the analysis of the free-free beam and coupled center load reso- 
nant vibration syste^ it was assumed that the exciter moving element^ 
the test cylinder and its supporting yoke and the hardware used in coup- 
ling the exciter moving element to the beam acted as a concentrated mass 
of magnitude m^at the beam center. The beam had a length and a 
total mass rrij^ . The coordinates used in the analysis are as shown below. 



The general equation for the transverse vibration of a prisimatical 



X = 




X = X 



‘n 



y 

y 



beam, as presented by Timoshenko 



(13) , 



is: 




(15) 



71 . 



which may be written: 



where 



E 

I 

§0 






dt^ Ar ^ 



0 



2 

modulus of elasticity, Ib/in^ 

beam cross-sectional moment of inertia, in 

2 

gravitational constant, 386 in/sec 

O 

beam weight density^ Ib/in*^ 



(16) 



Considering only the first mode of vibration as was justified in 
Chapter III, Section 1-b, the solution to Equation (l6) has the form 
y = X sinuOt , where sinoot gives the periodic variation of the 
beam deflection along its length according to the function X = f(x). 
Substituting this solution into Equation (l6) yields: 



d^X 



= 



(IT) 



where o( is defined by: 

OC** = ^ (18) 

Particular solutions of Equation (l8) involve the functions sint?')^ , 
cos(// , sinh(?(/. and cosho^./ . Therefore, the solution to Equation (l6) 
can be conveniently expressed as: 



y = X sin cdt = jci (cos(XX + coshiAV ) + Cg (cospty. - coshoC/. ) 

+ (sin<XX + sinhoty ) + (siny/, - sinhyx )jsinDoir 

(19) 



n 



cffi;;? 






I 
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Application of the end boundary conditions, ice., the bending 

d^y (i3y ^ 

moment El — - and shear force El — ^ vanish at x = 0 and x = > 



dx^ 



dx3 



reveals that C2 = C|^ = 0. Similarly, the mid-span boundary condition 



of zero slope. 



dx 



= 0 at X 



-I . 



gives : 



cosOf'Y + coshoiJ^ 
sin<y/ - sinh^</ 



( 20 ) 



An additional condition at beam center is obtained from the force- 
inertia balance at x = yf : 



(a) Exciter force (up): 

(b) Shear force (up): 

(c) Inertia load of m^ (down): 



F = 



Fq sinoo"t 



2EI 



d^y 

dx3 



- m 



<i^y 

dt2 



Thus, 



at X ^ J : 



F 

o 



sin 



U)~t 



+ 2EI 





(21) 



Substitution of Equations (19) and (20) into Equation (2l) results 
in the following: 



F + 2EI (X C3 



- cos 



C./f + cosho^/ + (sin + sinhoC>^ ) 



+ coshCl-f 



sin«/ - sinh 
{cos^i + cosho(./ )^ 



= - ^ 



I 0 n >.cosot^ + cosno..! 

S L sln«(^ - sinh«^ 



( 22 ) 









I 



I 
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Also, from the definition of the beam mass, nij^, may be conveniently 
introduced into Equation (22) as follows: 



, 2EI (ATw!) = 

Elg^ 2 </J 



( 23 ) 



For convenience, the first and second terms in brackets of Equation (22) 
are called Q and R. Substituting Equation (23 into (22) gives: 



- F 



M) 

c 



Q ( 









(24) 



Examination of Equations (l9) and (22) shows that beam displacement, 
y, is a function of and Of/ . At the theoretical resonance condition^ 
y becomes infinite. However^ for a finite frequency and a beam of finite 
mass and dimensions, 0(J( cannot approach infinity. Hence, with regard to 
theory, must become infinite at resonance. Therefore, the denominator 
of the right-hand side of Equation (24) must vanish at resonance, and; 





( 25 ) 



Substitution of the bracket terms Q and R and simplifying produces the 
following frequency equation in terms of beam to center mass ratio at 
resonance: 



. ^ (sech od + cos och 

m (sin + COS oU( tanh <U) 




74 . 



Equation (26) is plotted for various mass ratios versus values of 0(^ in 
Fig. 6. It is worthwhile to note that the first mode solution of the 
simple free-free beam with no added center-span mass should correspond to 
the value of (yj as 
frequency solution for this case as: 



m^ approaches infinity. Reference l4 gives the 



Elg 



tO = 22.4 



(27) 



AY'(2i’)^ 

Substitution for d into Equation (27) and solving for 0<-/ gives a 
value of oJ = 2.37 at resonance. Fig. 6 clearly indicates that 
approaches 2.37 as approaches infinity. 

In order to have minimum force transmission from the beam to the 
beam supports^ these supports must be located at the vibration npdal 
points. Solving for beam deflection, y, by substituting Equation (20) 
into ( 19 ) : 



y = 



cosO(i* + coshni^ , , X , , , ,T1 

- sinhi*/ 



sino)t = 0 



(28) 



At the nodes, say at x = x^, y = 0. Equating Equation (28) to zero gives 



sin 



od 



X/1 

J 



+ sinh(X/ — 



cos 



+ cosh C(~f 



cos 






+ cosh o(A 

A 



sinh^/ - sin 



(29) 
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The ratio of nodal distance from beam end to beam half length, ■> 

plotted versus in Fig. 7. 

The maximvun bending moment in the beam will occur at x = when 

the center amplitude, a, is at its maximum value. This will take place 

at times when sin is unity. From elementary beam theory, the bending 

moment is given by M = El — ^ ; thus, differentiating Equation (28) twice: 

dx 



M = El — ^ 
max dx 






2 sin sinhX^f , 

= 2E10( ^sln^yf -sinh^^~) (30) 



The maximum stress, ^max > given by times h/2 where h is the 

beam height in the direction of deflection. Combining this with Equations 
(24 and ( 28 ) gives: 



F h 1 
max " 2m 



r 



u. 



^ 



(31) 



where 



U, 



U, 



U 



= sin lyj tanh (32) 

= sin C/J + cos tanh Oi-f (33) 

^ = ^./(sechfl^/ + cos 0^.^ ) (34) 

Ug and are plotted versus in Fig. 8. 

The maximum amplitude at x is given by Equation (28) for times 

when sin d is unity. Upon substitution for C^: 
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F 1 

a = — —T 2 



U 



m CC>^ m 

=^"2 ^ “3 

c 



(35) 



Equations ( 31 ) and (35) are not useful at resonance because the 

m, 

denominator term — + U vanishes at the resonant condition. This 

m 3 

c 

is a result of the earlier assumption that maximum amplitude goes to 
infinity at resonance, i.e., approaches infinity at resonance. However, 
the ratio of bending moment to amplitude may be obtained by dividing 
Equation ( 3 I) t»y (35)^ and; 



cr 



max 



E _h^ 

2 



/ 



U, 



u. 



(36) 



The maximum bending stress versus amplitude and vibration intensity, af, 
for steel and aluminum beams at various a/ values and at f = 100 cps is 
plotted in Fig. 9* 

The interpretation of this analysis and final beam selection is 
presented in Chapter Section 1-b. The curves developed by means of 

this analysis (Figs. 6^ 8 and 9) were used for final beam design ac- 

cording to the calculation method indicated in Appendix B. It is con- 
sidered that this analysis served as a remarkably useful tool in this 
project. The resonant frequency obtained of the actual beam was vithin 
1 1/2 per cent of the design value. 
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APPENDIX B 

RESONANT BEAM DESIGN CALCULATIONS 



The resonant free -free beam design calculations based on the 
analysis presented in Appendix A are given below. Beam wi(ith, b, was 
initially set at I.50" so that existing hardware for test cylinder sup- 
port could be employed. The material selected, "HY- 100 " steel, had a 
weight density of O.283 Ib/in^. As explained in Chapter II, Section 1-b, 
a value of (/J = 2.25 was selected for the final design on the basis of 
stress and beam size considerations. The weight of the center span load, 
consisting of the test cylinder and test cylinder yoke, exciter extension 
and coupling hardware, and the exciter moving element, was estimated to 
be = 5*5 Its . Thus , for 0 (J ‘ 2.25 

“b ^b 

— - ^ = 5 -T 6 from Fig. 6 

and 



= 5.76 X 5.5 = 31.68 lbs 



also 



= b X h X 2 
b 



/ xy. 31.68 lbs 

Another equation foryf in terms of h may be obtained from Equation (18), 



Appendix A: 



2/2 
r AroO 12 V , 2TlT 

0 ( = ElgQ ” Eg^ V h ) 



Then 



0 (y( = 



^ (2JlfN2 

Eg,_ h 



1/4 

/ 



2.25 



I 






Squaring both sides of this equation gives: 



^ ^ 1/2 / 

<i^) (2Hf) 4- 



= (2.25) 



where 



f = 
E = 



So = 



100 cps 

6 . 2 

30 X 10 Ib/in 

2 

386 in/sec 



thus 



,12 X 0 



.283,^/^ 



( Z— ) (2TT X 100) 

30 X 10° 



L 



h 



(2.25) 



4 " = ^^5, or h = ^5 

Substitution of this relation for h into the above expression for W 



gives : 



W, = b X h X 2 
b 



/ X r = 31 >68 = b X 

475 



/ 



Thus may be calculated: 



(j 

1.5 X X 2y{ X 0.283 = 31*68 



/ 

J 



Beam length 



-•>31.68 xJ^T 5 

1.5 X 2 X 0.283 

26.07" 

2 y = 52.14" 
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Since h 



J ^AT5: 

2 

(26.07) 

n = — ; 

475 



1.43" 



The steel beam dimensions for a resonant frequency of 100 cps with a 
center load of 5*5 ll^s are: 



1.50” X 1.43” X 52.14” 



The steel beam as machined had the following dimensions; 



1.4T5” X 1.408" X 52.20” . 



The resonant frequency of this beam was determined as follows: 






W 

c 

ni„ 



0.283 X 1.4 t 4 X 1.408 X 52.20 = 30.717 lbs 

5.50 lbs 



w 

_ 30-717 
■ 5.50 



5.585 



From Fig. 6 , for 




= 5-585, 0 (/= 2 . 247 ; 



hence : 







2.247 

26.10 



12 y 



h 



0.0861 

1/4 



0.0861 



80. 



12 X 0.283 
30 X 10^ X 386 



,2XVf 2 

^T7IJU8) 



= (0.0861) 



and 



f = 97-0 cps 



It was desired to have a minimum resonant frequency of 100 cps, 
however. The calculated resonant frequency was raised to IO5.3 cps hy 
shortening the team to a length of 5O.OO". Thus, the above calculation 
procedure is used again, as; 

w = 0.283 X 1.4 t 4 X 1.408 X 50.00 = 29.367 lbs 
b 

29 

5 



from which 

f = 105.3 cps. 



.367 _ 



243 

.00 



5.340, = 2.243 from Fig 6 

1/4 

.0897 = ot ^^‘,2 






h 



The actual experimental resonant frequency of the system was 104 cps . 

The nodal distance from beam end, x^, was determined from Fig. 7 
with the value of Q(yf = 2.243: 

for ^ 1 ^ 7 = 2.243, = 0.503 

y 

since ^ = 2 5. 00" 

X = 25 X 0.503 = 12.575" 



Each gimbal assembly was then placed this distance from each beam end. 







I 



I 
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APPENDIX C 

CALIBRATION PHASE HEAT TRANSFER SAMPLE CALCULATIONS 

The calculations presented below were those performed for calibra- 
tion phase Test No. 9 (Table l) which had the following relevant original 
data: 

Q = watts 

m 

t = 266. 7°F 
s 

t = 76 .T°F 

a 

a. temperature potential 
Formula; At = t - t 

S £t 

t = 266. 7°F 
s 

t = 76 . 7 °F 

a 

At = 266.7 - 76.7 = 190. 0°F 

b. 0^, power delivered to test section 

Formula: = O.989I 0 ^ 

= 45.00 watts 

= 0.9891 X 45.00 = 44.510 watts 

c. 0^, heat transferred by radiation 

Method: = value read from Fig. 12 , corresponding to A t 

A t = 190. 0°F 

= 0.433 watts 



I 

I 

I 
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d. Q , heat transferred by free convection in absence of vibration 

CO 

5 A 

Formula: Q = 0.01428 ^ t , watts 

CO 

t = 190. 0°F, = 686.804 

Q = 0.01428 X 686.804 = 9.809 watts 

CO 

e . ^ heat loss from test cylinder 

Formula: ^o = ^ ’ ^^r %o^ 

0 ^ = 44.510 watts 

Q = 0.433 watts 

r 

Q = 9*809 watts 

CO 

Qd + Aio = ^^*510 - (0.433 + 9.809) 

Qd + Q^q = 3^*268 watts 



}* 



I 
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APPENDIX D 

VIBRATION PHASE HEAT TRANSFER SAMPLE CALCULATIONS 

The calculations presented below were those performed for vibration 
phase Test No. 6 (Table 2) which had the following relevant original data: 



f 


= 104 cps 


Q 


= 44.90 watts 


m 


a 


= 0.0836" 


^ ^hm 


= 52.11°F 


^s 


= 228.03°F 


t 

a 


= 78.37°F 



a. temperature potential 

Formula: At = t - t 

s a 

tg = 228.03°F 

t = t8.3T°f 

a 

= 228.03 - 78.37 - i 49.66 °f 

power delivered to test section 
Formula: = O.989I 0^ 

0^ = 44.90 watts 

= 0.9891 X 44.90 = 44. 4l watts 

c. heat transferred by radiation 

Method: = value read from Fig. 12, corresponding to t 

At = i49.66°f 

= 0.314 watts 




4 
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d. Qqo ^ heat loss from test cylinder 

Method: = value read from Fig. 13^ corresponding to 



/v t = 52.11°F 

hm 

+ Q^o "" 31-63 watts 

e . y heat transferred by convection in presence of vibration 

Formula: . (q^ + 

0^ = 44.i|-l watts 

= 0.31^ watts 

^ 31*63 watts 

= kh.kl - 0.31 - 31.63 = 12.47 watts 

f. h^^ heat-transfer coefficient in presence of vibration 

Q 

_ . 

Formula: h = 3.412 T . IT 

V A t 

Q = 12.47 watts 

cv 

A = 0.0993 

^ t = i 49.66°F 



12.47 




h 



0.0993 X 149.66 



V 



g. af, the intensity of vibration 



Formula: af = 

a = 0.0836" 



af 



f 



.104 .cps 
0.0836 X 104 



0.725 ft /sec 
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